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Abstract 

Minimum Bias production cross sections of g mesons have been measured in 200 AGeV/c 
S+Au and S+S collisions at the CERN SPS by reconstructing the 77 — 7-77 decay. 
The measurements have been made over the rapidity range 2.1 < y < 2.9 using the 
leadglass spectrometer of WA80. Within the statistical and systematical uncertain¬ 
ties the spectral shapes of 7 r° and 77 mesons yields are identical when their invariant 
differential cross section is plotted as a function of the transverse mass. The relative 
normalization of the 77 to 7 r° transverse mass spectra is found to be 0.53 ± 0.07 for 
S-fAu and 0.43±0.15 for S-fS reactions. Extrapolation to full phase space leads to an 
integrated cross section ratio of 77 to 7 r° mesons of 0.147± 0.017(stat.) ± 0.015(syst.), 
and 0.120 ± 0.034(stat.) ± 0 . 022 (syst.) for S-fAu and S-fS collisions, respectively. 
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1 Introduction 


The measurement of inclusive transverse momentum distributions of produced particles 
and their relative production cross sections is known to be a valuable tool in the study of 
high energy hadronic interactions. While many efforts have been undertaken to measure 
charged hadrons in ultra-relativistic nuclear collisions at the CERN-SPS and BNL-AGS 
[1], up to now only little is known about the production of neutral mesons heavier than the 
7r°. The WA80 collaboration at the CERN-SPS has therefore concentrated on identifying 
7r° and 1 ] mesons with high precision over an extended transverse momentum range via 
their two-photon decay channel. Being the lightest mesons, with the t] meson containing 
(hidden) strangeness, they provide important information about the dynamics and decay 
of the hot hreball formed in central nucleus-nucleus collisions. 

In pp collisions, relations between the production cross sections for different mesons have 
been discussed within the phenomenological concept of scaling with transverse mass niT = 
ttIq p^. In this picture, the shapes of the differential cross sections of mesons, plotted 
as a function of niT are expected to be identical. As a simple consequence of different 
rest masses, this leads to a strong px-dependence of the cross section ratio at low values 
of pt (see e.g. [2, 3] and Pig. 5). Por values of px on the other hand, the meson 

yield ratios were found to become constant and independent of the reaction system and 
CM energy. Some qualitative interpretations of this behavior are discussed in the literature 
[4, 5, 6]. 

Previous data on mid-rapidity pion and kaon production in nuclear reactions at CERN [7, 8] 
and BNL [9] seem to be consistent with this picture. However, the dynamic range in these 
experiments generally spans only a limited transverse momentum range at low px and their 
statistical and systematic errors still leave some freedom. In the present work we will extend 
these studies to neutral pious and ?y-mesons. A possible scenario, in which a deviation from 
niT scaling may be observed, is an expanding hreball including collective transverse how, 
because in this picture thermal and how components ahect mesons with diherent masses 
diherently [10]. Therefore, in the presence of a large collective how component and similar 
freeze-out conditions, the mx spectra of heavier mesons are expected to become hatter. 
The present absolute cross section measurements of 7r° and t] mesons with approximately 
3 GeV !c overlap in transverse momentum and with identical rapidity acceptance provide 
a stringent test of the mj scaling concept. 

2 Experimental 

The WA80 experiment at the CERN-SPS (Eig. 1) is equipped with a hnely segmented lead- 
glass spectrometer covering the polar range 6.2° < diab A 13.9°. Compared to preceding 
WA80 data taking runs, the leadglass calorimeter SAPHIR [11] has been augmented by two 
additional leadglass calorimeters (TOWERS) in order to achieve an improved and uniform 
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acceptance for 7 r° and t] mesons in the rapidity range 2.1 < y < 2.9. Both meson types 
are reconstructed by invariant mass analysis of their 77 decay mode. The total azimuthal 
coverage of the 3798 leadglass modules at a distance of 9 m to the target is ~ 1.3 tt. The 
SAPHIR modules consist of SF5 leadglass blocks of dimension 3.5 X 3.5 X 46 cm^, whereas 
the TOWER modules are made of 4.0 X 4.0 X 40 cm® TFl leadglass. 

The minimum bias trigger requires a valid signal of the beam counters and a minimum 
amount of transverse energy Et ^ 1 GeV, detected by the Mid-Rapidity Calorimeter [12]. 
Data have been taken with the 200 AGeV/c sulfur beam on targets of An (250 mg/cm®) 
and S (205 and 510 mg/cm®). For the present analysis 7.9 million S-|-Au and 1.4 million 
S-|-S minimum bias events were accumulated. The px distributions for 7 ’s and 7 r°’s from 
S-|-Au interactions could be determined over a wide range of 0.3 < px A 4.5 GeV/c. For 
the hrst time also 7 meson yields are measured in ultrarelativistic heavy ion reactions over 
the pt range of 0.5 < px A 3.5 GeV/c. 


3 Data Analysis 

In addition to modihcations of the experimental setup more rehned analysis methods have 
been developed to improve the particle identihcation performance compared to previous 
WA80 analyses [13, 14]. Most importantly, a more efficient distinction between hadronic 
and electromagnetic showers is achieved by analyzing their lateral dispersion. Furthermore, 
an improved procedure for unfolding overlapping showers was used [15]. Different than from 
other works [16], the charged particle veto was not used in the present analysis. 

In order to reconstruct 7r° and 7 mesons from the experimental data and to extract their 
differential cross sections, dcr/dpx, the invariant mass of each 77 -pair of an event is cal¬ 
culated and sorted as a function of the transverse momentum of the photon-pair. The 
resulting invariant mass distributions are composed of peaks centered at values of the 7 r° 
and 7 rest mass and are located on top of a large combinatorial background. The meson 
yields are extracted from these distributions by subtracting the combinatorial background. 
Possible contaminations in the photon sample, e.g. due to non-rejected charged hadrons 
(see above), thus do only contribute to the uncorrelated (and subtracted) combinatorial 
background, but not contribute to the isolated invariant mass peaks itself. Due to the 
high particle multiplicity in heavy-ion collisions the resulting signal-to-background ratios 
may become very unfavorable particularly at transverse momenta below 1 GeV/c, so that 
precise knowledge of the shape and yield of the combinatorial background below the mass- 
peaks is required. Instead of htting the observed background distributions by arbitrarily 
chosen continuous functions, we have therefore applied an event-mixing procedure to cal¬ 
culate the background distributions from the experimental data itself [17, 18]. Special care 
has been taken to ensure that the data and background distributions are generated for the 
same class of events with very similar phase space characteristics. Therefore, photons from 
a given event class have been combined exclusively with photons from another event of 
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the same centrality class, which have been selected using 8 equidistant bins of the associ¬ 
ated transverse energy Et- To minimize the statistical and systematic uncertainties in the 
extraction of the t] peak contents, a two-dimensional ht of the data/background distribu¬ 
tion as a function of px and invariant mass has been performed. Based on the measured 
response function of the calorimeter, the peak has been assumed to be of Gaussian shape 
with position and width taken as free parameters but hxed to be the same for all px bins. 

To validate the procedure of backgound generation and to check the accuracy of the me¬ 
son reconstruction procedure, a Monte Carlo simulation was performed in which a known 
number of 7r° and 1 ] mesons were generated according to the experimentally observed dis¬ 
tributions. After subjecting those photons found within the experimental acceptance to 
realistic detector resolution effects, the same analysis techniques as used for the experi¬ 
mental data were applied. The number of mesons extracted agrees with the known number 
of generated particles [19, 18] within the statistical errors of the background subtraction 
procedure. Using similar statistics as in the experimental data, systematic deviations im¬ 
posed by the background subtraction were found to be about 2% for the 7r° and 10% for t] 
distributions. 

Fig. 2 shows 77 invariant mass distributions for the S-|-Au data in the 7r° and t] mass range 
after background subtraction. Whereas the signihcance of the 7r° peak is obvious, that of 
the 1 ] peak needs to be proven. All invariant mass distributions were thus htted in the 
range of the t] peak both by assuming a constant (ht a) and the sum of a Gaussian plus a 
constant (ht b). The ratio of the obtained chi-squares, r = xla/xlh^ yields the signihcance 
of the peak via an U-test [20]. In the case of the t] peak of Fig. 2b, the signihcance level is 
approx. 95 % and even the hts obtained for narrower px ranges result in signihcance levels 
between 70% and 90%. The position and width of the experimental peaks were found 
to be compatible with the known t] mass of = 547.45 MeV/c^ and with the expected 
invariant mass resolution, = 25 MeV/c^, taking into account the energy and position 
resolution of the leadglass calorimeters [21]. 

The 7r° and 1 ] distributions are corrected for the geometrical acceptance, which considers 
the loss of photons due to the incomplete coverage of the full solid angle and the software 
imposed low energy photon cut-oh of E^in = 750 MeV. Furthermore, the reconstruction 
efhciency has been determined in order to correct for losses during analyis of photons, and 
thereby 7r° and t] mesons, due to the shower reconstruction procedure. Such a correction is 
of particular importance in high multiplicity heavy-ion reactions where a large probability 
for overlapping showers is encountered which may alter their reconstructed energies, po¬ 
sitions, and transverse dispersions. To calculate the efficiency, GFANT simulated photon 
showers, which were well adjusted to the measured detector response, were superimposed 
onto real data events. By analyzing such artihcial events with the same package of shower 
reconstruction programs as used for real data, the reconstruction efficiency for those super¬ 
imposed particles could be extracted. A detailed description of this method [21, 19] will be 
presented in a forthcoming publication. Finally, we point out that effects of reconstruction 
efficiencies and other systematic errors largely cancel in the Ty/vr^-ratio as a function of my. 
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because both mesons were identified via the same 77 -decay mode and have similar energy 
for fixed transverse mass of the parent particle. The total systematic uncertainty of up to 
approx. 11 % in the extracted ^/vr® ratios thereby are dominated largely by the extraction 
of the 1 ] peak contents from the invariant mass spectra. 


4 Results and Discussion 

Acceptance and efficiency corrected invariant cross sections of 7 r° and 7 mesons are shown in 
Fig. 3 for minimum bias S-|-Au and S-|-S reactions as a function of the transverse mass niT- 
In addition to the statistical errors, the systematic errors of the acceptance and efficiency 
correction have been added quadratically. The invariant cross section of 7 r° mesons in both 
systems clearly exhibit a spectrum with concave shape. A similar shape is observed for the 
7 mesons. This so called niT scaling behavior of the invariant cross sections is supported by 
a variety of experimental results and has been observed for a wide range of center of mass 
energies and for different meson masses [7, 9, 22, 23, 24]. Phenomenological interpretations 
are mostly based on thermal particle production mechanisms and are discussed for example 
in Refs. [4, 5, 6 ]. 

In order to analyze the data in more detail, an empirical power-law formula (Eq. 1) proposed 
in Ref. [5] 

£ A = c (-^) 

\Po + J 

has been fitted to the 7 r° niT spectra of Fig. 3 with the free parameters (7, po, and n. For 
the S-|-Au data the extracted values are po = 6.8 ± 0.1 GeV/c and n = 37.5 ± 0.3. The 
results are displayed in Fig. 3 as solid lines. The dotted lines, which represent fits to the 
7 distributions assuming the same values for po and n as for the 7 r° fits, provide a good 
description of the data (y^/d.f. = 10.1 /5). A separate fit with 3 free parameters (dashed 
line) is almost identical in the range of the 7 data (y^/d.f. = 9.7 / 3). To further check the 
sensitivity of the data to niT scaling properties, we have analyzed the 7 and 7 r° invariant 
cross sections by fitting pure exponentials in the same px and mx ranges, respectively. It 
turns out that in the mx representation the inverse slope parameter of the 7 distribution, 
Tr],mT = (226 ± 9) MeV, is only 4 % (= Icr) larger, while in the px representation it is 15 % 
(^ 3 (t) larger than the corresponding values deduced from the 7 r° spectra. These numbers 
demonstrate the sensitivity of the data and seem to indicate at most a very small deviation 
from mj scaling. In Ref. [10] the energy distribution of particles is interpreted as being 
composed of a thermal and a collective radial expansion component. Applying this picture 
to the observed (but not significant) deviation of the 7 r° and 7 mj slopes suggests a radial 
expansion velocity at freeze-out of appr. 30-50 % of the speed of light, well in line with the 
results of Ref. [10]. 

To test the concept of mx-scaling in a model independent way, the differential 7 / 7 r° ratio is 
presented in Fig. 4 for identical bins in mx. The average of all data points yields a ratio of 


n 

( 1 ) 
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0.53ib0.07 for the S+Au data and 0.43±0.15 for the S+S data. Within the error bars and 
over a broad niT range we observe no signihcant deviation from a horizontal line, again 
demonstrating consistency with niT scaling. 

Figure 5 shows as an alternative representation of the measured data the ?y/7r°-ratios as 
a function of px- The full symbols represent the results obtained within the present and 
previous experiment [14] for three different reaction systems (see also Table 1), whereas the 
open symbols show a compilation of various p+p, vr+p, and p+C midrapidity [25, 26, 27, 
28, 29, 31] and forward rapidity data [2, 3, 30]. The latter two have predominantly been 
taken at higher CM-energies (-^ = 19.4 - 62.0 GeV) as compared to the present experiment 
(-^ = 19.4 AGeV). A comparison of all data points might nevertheless be justihed, since 
the ?y/7r°-ratios have been found to be insensitive to ^/s in p-induced reactions [26]. For 
a quantitative comparison, the spectral shapes of the different data sets were taken into 
account again according to Fq. 1. Plotted as a function of px the differential ?y/7r°-ratio is 
then written as 


A — 
dp^ 


id) 


^d^a , r 


= R 


'Tj/ti 



+ Px 


+ Pt, 


( 2 ) 


with po and n as in Fq. 1 and Rni-K = Cnj Using po and n as extracted from the 7r° mx 
distributions, the differential t] j cross section ratio is found to be = 0.55 ± 0.07 

in S+Au data (solid line in Fig. 5) and = 0.42 ± 0.13 in S+S data, well in agreement 

to the results above. For comparison, the differential ?y/7r°-ratio of the available p-induced 
data was extracted in the same way. Taking as an average po = 4.9 ± 0.1 GeV/c and 
n = 33.6 ± 0.7 from pion data of various experiments [25, 32, 33, 34, 35, 36] we hnd 
= 0.55 ± 0.02 (dashed line in Fig. 5). 

It should be emphasized, that the numbers quoted above represent the cross section ratios 
measured for equal bins in nix- Because mx ~ px for px rriQ^ the same result would be 
obtained for large values of px- Fxtracting the integrated production ratio requires 

an extrapolation of the experimental px spectra to px 0 , or likewise an extrapolation 
of the experimental mx spectra to mx —> m,f^ and m^, respectively (see Fq. 2 and Fig. 5). 
Such an extrapolation into regions outside the experimental acceptance induces additional 
systematic uncertainties. Assuming the ht function of Fq. 1 the extrapolation yields a ratio 
of integrated production cross sections of [g j tt^)s+Au = 0.147 ± 0.017(stat.) ± 0.015(syst.), 
(?7/7r°)5+5 = 0.120 ± 0.034(stat.) ± 0.022(syst.), and (?7/7r°)p = 0.103 ± 0.003(stat.) ± 
0.006(syst.) for the different reaction systems, respectively. As a consequence of the larger 
inverse slope of the spectra in S induced reactions, the integrated gjir ratios appear to be 
enhanced over p induced reactions, even though the asymptotic ratios Rni-K are very similar. 
For comparison with other mesons, the K'^ ratio in central 14.5 AGeV Si+Au collisions 
is measured to be 0.185 ± 0.007 [9], while a value of 0.137 ± 0.008 has been deduced from 
central 200 AGeV/c S+S collisions [37]. Similarly, A"°/7r“ = 0.104 + 0.012 may be inferred 
from 200 AGeV/c S+Ag data [8]. Fmploying a simplihed quantum statistical hadron- 
gas model in which relative particle production of sufficiently large systems is governed 
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reaction 

system 

px-range 

(GeV/c) 

O 

o' 

mx-range 

(GeV/c) 

O 

o' 

S+Au 

0.5 - 0.7 

0.175 ±0.066 

0.7 - 0.9 

0.512 ±0.162 


0.7 - 1.0 

0.163 ±0.065 

0.9 - 1.1 

0.250 ±0.149 


1.0 - 1.5 

0.423 ±0.087 

1.1 - 1.5 

0.674 ±0.147 


1.5 - 2.0 

0.398 ±0.115 

1.5 - 2.0 

0.461 ±0.148 


2.0 - 2.5 

0.705 ±0.183 

2.0 - 2.5 

1.048 ±0.232 


2.5 - 3.5 

0.496 ±0.203 

2.5 - 3.5 

0.581 ±0.231 

S+S 

0.5 - 1.0 

0.219 ±0.074 

0.7 - 1.1 

0.510 ±0.194 


1.0 - 1.5 

0.173 ± 0.122 

1.1 - 1.5 

0.306 ±0.235 


Table 1: Measured ?7/7r° cross section ratios for different bins in px and niT- 

by their respective masses and in which modihcations due to resonance decays are taken 
into account, a temperature larger than T ^ 150 MeV is needed to reproduce the present 
experimental data [38]. 

In summary, inclusive differential cross-sections of t] mesons have been measured for the hrst 
time in 200 AGeV/c nuclear collisions. Together with the 7r° spectra measured within the 
same apparatus via the same decay mode, the results are compatible with the phenomeno¬ 
logical concept of niT scaling. The slight deviation observed, can be well accounted for by 
assuming a radial expansion velocity of /3 sa 0.3-0.5. The ratio of normalizations of the niT 
spectra, or equivalently, the ?y/7r°-ratio at large px, were found to be 0.53 ± 0.07 for S-|-Au 
and 0.43 ± 0.15 for S-|-S collisions. Extrapolation to full phase space yields an integrated 
production cross section ratio of 0.147±0.017(stat.)±0.015(syst.), and 0.120±0.034(stat.)± 
0.022(syst.), respectively, which is to be compared to 0.103 ± 0.003(stat.) ± 0.006(syst.) 
extracted from data of p-induced reactions. 
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Figure 2: Background subtracted invariant mass distributions 
mass range accumulated for transverse momenta of 0.5 < pr ^ 
background ratio is 5.7 % for (a) and 0.07 % for (b). 
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Figure 3: Invariant cross sections of 7r° and t] mesons averaged over the rapidity interval of 
2.1 A J/ A 2.9 as a function of transverse mass for 200 AGeV/c S+Au and S+S minimum 
bias data. The S+S data have been scaled by a factor 0.5 for better presentation. The 
different lines are htted results as explained in the text. 
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Figure 4: ?y/7r°-ratio as a function of niT- The lines represent the mean values of the S+Au 
(solid line) and S+S data points (dashed line). 





Figure 5: ?y/7r°-ratio as a function of px- The solid and dashed lines represent nix scaling 
parameterizations (Eq. 2) extracted from S+Au data and p induced collisions, respectively. 
(See text for details.) 
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